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Abstract. 

While multi-layer systems can possess global inversion centers, they can have regions with 
locally broken inversion symmetry. This can modify the superconducting properties of such 
a system. Here we analyze two dimensional multi-layer systems yielding spatially modulated 
antisymmetric spin-orbit coupling (ASOC) and discuss superconductivity with mixed parity 
order parameters. In particular, the influence of ASOC on the spin susceptibility is investigated 
at zero temperature. For weak inter-layer coupling we find an enhanced spin susceptibility 
induced by ASOC, which hints the potential importance of this aspect for superconducting 
phase in specially structured superlattices. 



1. Introduction 

The discovery of superconductivity in the non-centrosymmetric heavy Fermion compound 
CePtsSi [U [2] has attracted much attention. Non-centrosymmetricity leads to antisymmetric 
spin-orbit couphng (ASOC) which, for example, gives rise to mixed parity pairing, and the 
characteristic anisotropy of the spin susceptibility in the superconducting state [3-9]. 

Recently, artificial superlattices involving the heavy Fermion compound Celns [lOJ and 
CeCoIns [llj have been fabricated, realizing two-dimensional multi-layer structures of heavy 
electron systems and ordinary metals. A particularly interesting case is the superlattice of 
CeCoIns and YbCoIns which shows superconductivity, most likely induced by the CeCoIns 
layers. Motivated by this system we setup a model of a superconductor with similar layered 
structure, including the fact that such structures show distinctive violation of inversion 
(reflection) symmetry, while the overall system has inversion centers. Such local non- 
centrosymmetricity is expected to have a pronounced effect on the superconducting phase 
through local occurrence of ASOC, in particular, for the spin susceptibility. 

2. Formulation of model 

We introduce here model Hamiltonian for multi-layer systems including ASOC, given by 
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where c^sm ((^ksm) annihilation (creation) operator for an electron with spin s on the 

layer m, and a^s' is the Pauli matrices. The (x, y, 2;)-axes correspond to the (a, b, c)-axes of the 
tetragonal crystal. 

The first two term describe the intra-layer dispersion and the ASOC. We consider a square 
lattice with a tight-binding model, i.e., £(k) = —2t{coskx + cos ky) — jj.. We choose the 
unit of energy as t = 1 and assume the chemical potential ^ = —1. The electron density 
per site is approximately 0.63. The ASOC term preserves time reversal symmetry through 
the condition g(— k) = — g(k) and has Rashba structure, which we assume a simple form 
g(k) = {— sin ky, sin kx,0). The coupling constants a„j arc layer dependent and have opposite 
sign for layers above and below a center layer. The third term introduces intra-layer Cooper 
pairing via an off-diagonal mean field. Here Ass'mi^) involves both the spin singlet and triplet 
components, 
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where V'm(k) and dm(k) are scalar and vector order parameters for the spin-singlet and triplet 
pairing on layer m, respectively. For simplicity we use an order parameter on phenomenological 
grounds without resorting to any microscopic model based on a pairing mechanism. To be 
concrete we use an s-wave order parameter for the singlet and a p-wave order parameter for 
the triplet pairing, where on symmetry grounds we request df„(k) || g(k): Vm(k) = '0m and 
dm(k) = (im,g(k) = d„i{~ sin ky , sin kx , 0) ■ Wc choose \ijjm\, \dm\ < 0.01, small enough to satisfy 
the condition |As3'^(k)| <C la^l <S £f- The dominant order parameter component keeps the 
same sign over all layers, while the other (subdominant) component adjusts the sign with the 
ASOC ( 

o-rn)- The fourth term describes the intcr-layer hopping of electrons between nearest- 
neighbor layers. We assume that the inter-layer hopping t± is smaller than the intra-layer 
hopping t. 

3. Numerical results 

We now calculate the spin susceptibility of the multi-layer superconductors with spatially 
inhomogeneous ASOC, concentrating on the magnetic field direction along the c-axis. The spin 
susceptibility x = ^^t^h^o{Ms) / H is obtained numerically by calculating the magnetization 
{Ms) for a small magnetic field H. The necessary Zeeman coupling term is given by, 

Hz = -^ ^-(^ss'cls^Cks'm, (3) 
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where g = 2 and fj,B is the Bohr magneton. First, the Hamiltonian is diagonalized in the presence 
of a field, introducing the unitary transformation = f j^C/^(k) in Nambu-space of M layers, 
where the quasiparticle operators form a 4M-dimensional vector 

and analogous for the Bogoliubov quasiparticle operators f = (7^15 7k2' ■ ■ • >7k4M)- Thus, the 
Hamiltonian is 
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where -Ej(k) are the quasiparticle energies. The magnetization is obtained as, 
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where f{E) is the Fermi-Dirac distribution function. The matrix representation of spin operator 
is defined in the f basis as 

5^(k) = U\k)t''U{k), (7) 

with T,^ the /i-component of the spin operator in the 4M-dimensional space. 

As concrete examples we discuss the spin susceptibihty of bi-layer (M = 2) and tri-layer 
(M = 3) systems at T = 0. The corresponding couphng constants of ASOC are described as 
(ai,a2) = {a,— a) for bi-layers and (01,02,03) = (a, 0, — a) for tri-layers. 

We compare now the two cases: (1) the spin triplet channel is dominant \dm\ > \i^m\ and 
(2) the spin singlet channel is dominant |(im| < IV'ml- In case (1) the spin susceptibility remains 
unaffected by the superconducting state, Xs = Xn, because the spin triplet component of the 
type dm(k) oc g(k) _L z is an equal-spin pairing state with Cooper pair spins along the c-axis. 
Thus, spin polarization in the superconducting phase is possible without pair breaking. This 
feature is essentially independent of ASOC and inter-layer hopping as can be seen in Fig{T] for 
both the bi- and tri-layer systems. 




Figure 1. Spin susceptibility along c-axis for bi-layer system (i) and tri-layer system (ii). We 
assume t± = 0.1 (solid line) and t± = 0.2 (dashed line). 
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Figure 2. Spin susceptibility 
along c-axis for tri-layer system. 
The contribution from each layer is 
shown. We assume |V'm| > l^ml, 
and tx_ = 0.1 (solid line) and = 
0.2 (dashed line). 



More interesting is case (2) as spin singlet pairing leads to complete suppression of the spin 
susceptibility at T = in a conventional superconductor. Indeed for vanishing ASOC (a = 0) 
we find Xs = irrespective of t±. As soon as ASOC is turned on, however, the spin susceptibility 
gradually recovers and approaches for large a a constant value: Xs Xn for the bi-layer and 
Xs — ^ 2xn/3 for the tri-layer. The mechanism for this behavior lies in the spin-splitting of the 
electronic spectrum due to the Rashba-type ASOC. The spin polarization involves now an inter- 
band Van-Vleck-type contribution which is only weakly affected by the opening of quasiparticle 



gap in the superconducting phase. Note that this inter-band contribution rehes on the ASOC for 
matrix elements with the Zeeman couphng and is only available for the layers with non- vanishing 
a. Consequently, in the bi-layer system all layers are involved, giving rise to full recovery of 
Xs for large a (analogous to the uniformly non-centrosymmetric superconductor [5j), while in 
the tri-layer system only two of three layers can contribute yielding a correspondingly reduced 
limiting value for Xs- Figure [2] corroborates this picture by considering the contributions of the 
different layers. Indeed in the large a regime the outer layers m = 1, 3 carrying ASOC saturate 
at Xs Xn/3 while the center layer m = 2 completely suppresses. Remarkably at small a 
(< t±)i Xs behaves for all layers in the same way and leads for the center layer to a striking 
non-monotonic a-dependence. 

The numerical data in FiglHshow that the inter-layer hopping is in competition with ASOC, 
such that a larger t± yields a higher effective a^s ~ t± for the crossover from the behavior of 
conventional superconductor to that of non-centrosymmetric superconductor. This crossover is 
best evident in the peak of Xs around a^s ~ t± for the center layer (Figl2|). Thus, modifying 
G-g-) by uniaxial stress along the c-axis, can influence the magnetic response for c-axis fields 
in case (2). No such effect is expected for case (1). 

Within our model we find that the spin susceptibility along the a6-axis is always the half 
of the value observed along c-axis, independent of the strength of a and t± and the number 
of layers. Furthermore, we find that the spin susceptibility for both field directions is affected 
by the phase difference of order parameter between layers, but independent of the ratio of spin 
singlet and triplet components. Details will be explained elsewhere. 

4. Conclusion 

In view of the fact that CeCoIns is known to realize spin singlet superconductivity, we believe 
that most likely case (2) of our discussion is relevant for the multi-layer systems. Thus, the large 
observed upper critical fields in the superlattice of CeCoIns [llj would then rely on the presence 
of the spatially modulated ASOC. Moreover we believe that the variability of the superlattices 
and also the possibility of local measurements of magnetic properties through NMR would give 
many intriguing insights into the aspect of ASOC in these artificial systems. 
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